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Uncertainties in the Solar Neutrino Flux 
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I discuss three issues relevant to solar neutrino flux measurements: cross section uncertainties in pp chain 
reactions, uncertainties in the GALLEX/SAGE response to 7 Be and 51 Cr neutrinos, and the implications of 
helioseismology for nonstandard suns with mixed cores. A few comments are also offered on i/ e «-> v T oscillations, 
cosmologically interesting neutrino masses, and recent proposals for supernova neutrino observatories. 



1. INTRODUCTION 

It is a pleasure to be present for this historic 
meeting hosted by the Superkamiokande collabo- 
ration. In this talk I will address three issues af- 
fecting the solar neutrino flux and one connected 
with future detectors for supernova neutrinos. 

1.1. Nuclear Physics of the pp Chain 

One of the crucial inputs into the solar model 
is the network of nuclear reactions comprising the 
pp chain (and CNO cycle) . This network involves 
nonresonant charged particle reactions occurring 
at center-of-mass energies well below the height of 
the Coulomb barrier. As the solar core tempera- 
ture T c ~ 1.5-10 K, the typical kinetic energy for 
a nucleus in the core is < E >~ 2 keV. The com- 
petition between the Coulomb barrier and Boltz- 
man distribution leads to a typical energy for re- 
acting nuclei of (-^reacting) ~ 10 keV, a value that 
is generally lower than that where such reactions 
can be measured in the laboratory. Thus the task 
for nuclear physicists is to measure such reac- 
tions as accurately as possible over the accessi- 
ble range of laboratory energies, then extrapolate 
these measurements to the energies relevant for 
the sun. 

In the case of the driving reaction of the pp 
chain 



P + P 



Tl 



(1) 



the cross section is not measurable in the lab- 
oratory. Thus we must rely on theory. Fortu- 
nately deuterium is the simplest nucleus, and its 
properties are very well reproduced by NN poten- 



tials (Bonn, Paris, Argonne vl8, etc.) carefully fit 
to phase shifts. The calculated cross section de- 
pends on the accuracy with which the axial vec- 
tor coupling gA is known and on two-body cor- 
rections to the space-like component of the axial 
current, which are fortunately of order (v/c) 2 ~ 
1 %, where v is a typical bound nucleon velocity. 

The other major reactions of the pp chain are 
measureable, but generally not at the low energies 
relevant to our sun. The necessary extrapolation 
of the cross section a(E) to lower energies is ac- 
complished via the S factor 



a (E) = S ^ c -2*ZiZ**/P 
y ' E 



(2) 



where E is the center-of-mass energy, Z\ and 
are the charges of the interacting nuclei, and (3 
is the relative velocity. The introduction of S{E) 
removes the s-wave Coulomb interaction of point 
particles and thus provides a much smoother 
quantity for use in extrapolating data. S(E) de- 
pends on a number of physical effects - nuclear 
finite size, atomic screening corrections, higher 
particle waves, etc. - that the theorist must eval- 
uate before this extrapolation can be done. 

Because the solar neutrino problem is at a cru- 
cial juncture, a group of about 40 experts recently 
met at the Institute for Nuclear Theory, Seattle, 
to discuss the nuclear physics of the pp chain and 
CNO cycle. The questions addressed included the 
best current values for cross sections, critiques 
of assigned uncertainties, and recommendations 
for future experimental and theoretical work that 
could further improve our understanding of the 
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nuclear physics. The summary of this workshop 
will appear in Reviews of Modern Physics (Oc- 
tober, 1998) and is also available on the LANL 
preprint archive [1]. 
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Figure 1. The 7 Be(p, 7) 8 B S-factor as measured 
by Filippone et al. [3] and by Kavanagh et al. [4]. 
For each data set two theoretical extrapolations 
to S(0), reflecting different choices for the strong 
potential, are shown [2]. 



While I cannot give an adequate summary here, 
I will mention two of the reactions where signif- 
icant changes were recommended. The first of 
these is 7 Be(p, 7) 8 B, where the standard Si7(0)^ 
22.4 eVb is that given by Johnson et al [2]. 
Measurements of Si7(E) are complicated by the 
need to use radioactive targets and thus to deter- 
mine the areal density of the 7 Be target nuclei. 
Two techniques have been employed, measuring 
the rate of 478 keV photons from 7 Be decay or 
counting the daughter 7 Li nuclei via the reaction 
7 Li (d,p) 8 Li. The low-energy data sets [3,4] for 
Si7(E) disagree by 25%, a systematic effect ap- 
parent in Fig. 1 . Each data set is consistent with 
theory in its dependence on E: this dependence 
is simple in the illustrated low energy region as 
it is determined by the asymptotic nuclear wave 
function. 

The Seattle working group on Si7(E) found 
that only one low-energy data set, that of Fil- 



ippone et al. [3], was described in the published 
literature in sufficient detail to be evaluated. The 
target activity in that experiment had been mea- 
sured by both 478 keV gamma rays and by the 
(d,p) reaction, with consistent results. The re- 
sulting recommended value was thus based on this 
measurement, yielding 



517(0) = 191$eV b, la . 



(3) 



Since the workshop, two developments have 
occurred. The Orsay/Bordeaux/Paris-Sud/ 
USTHB group published [5] a new measurement 
of S17 (0) = 18.5 ± 1.0 eV b, while a preliminary 
value from the Weizmann/Troitzk/Mainz/Isolde 
group of S 17 (E p = 1.2 MeV) = 22.5 ± 2.5 eV b 
has been announced [6]. 

The 3 He(a,7) 7 Be reaction has been measured 
by two techniques, by counting the capture 7 
rays and by detecting the resulting 7 Be activity. 
While the two techniques have been used by sev- 
eral groups and have yielded separately consistent 
results, the capture 7 ray value Si7(0) = 0.507 ± 
0.016 keV b is not in good agreement with the 7 Be 
activity value 0.572 ± 0.026 keV-b. The Seattle 
working group concluded that the evidence for 
a systematic discrepancy of unknown origin was 
reasonably strong and recommended that stan- 
dard procedures be used in assigning a suitably 
expanded error. The recommended value S34 (0) 
is 0.53 ± 0.05. 

These and other recommended values were 
recently incorporated into the Bahcall and 
Pinnsoneault (BP98) solar model calculation [7]. 
While the workshop's recommended values in- 
volve no qualitative changes, there is some broad- 
ening of error bars and a downward shift in 
Si 7 (0), leading to the lower BP98 8 B flux. The 
workshop's Reviews of Modern Physics article 
summarizes a substantial amount of work on top- 
ics not discussed here: screening effects, weak ra- 
diative corrections to and exchange current effects 
on p+p, the atomic physics of 7 Be + e~, etc. 
Much of this discussion was useful in evaluating 
possible uncertainties in solar microphysics, and 
in identifying opportunities for reducing those un- 
certainities. 
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1.2. The Nuclear Physics of the 

G ALLEX / SAGE 51 Cr Calibrations 

The 51 Cr neutrino source experiments provide 
an important check on the overall gallium de- 
tector operations under few atom, hot chemistry 
conditions. The issue discussed here, and which 
was mentioned in the earlier experimental talks, 
is the potential complication due to contributions 
of uncertain strength to the 5/2~ and 3/2~ 71 Ge 
excited states (see Fig. 2). 
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Figure 2. Level scheme for 71 Ge showing the ex- 
cited states that contribute to absorption of pp, 
7 Be, 51 Cr, and 8 B neutrinos. 



The results of the source experiments can be 
normalized to the known 71 Ge ground state con- 
tribution [8], yielding 



Ro — E 



1 + 0.67 



BGT(5/2- 
BGT(gs) 



0.22 



BGT(3/2" 



BGT(gs) 
(4) 



_ f 0.98 ± 0.08, GALLEX [9] 
~ \ 1.00 ±0.13, SAGE [10] 

where E represents any departure of the efficiency 
from the value determined from tracers and used 



by the experimentalists in their analyses. The 
dependence of the results on the unknown BGT 
values is explicit. Clearly the conclusion E ~ 1 
requires an independent determination that the 
unknown BGT values are much smaller than the 
ground state value BGT(gs). 

It had been assumed that forward-angle (p,n) 
charge exchange measurements determine the un- 
known BGT values 



BGT(5/2~) 
BGT(gs) 

BGT(3/2~) 
BGT(gs) 



<0.06 



0.12 ±0.02 



(5) 



But extensive investigations [11,12] of the propor- 
tionality between (p,n) cross sections and known 
weak interaction BGT values have shown that 
the relationship is a complicated one. Existing 
discrepancies can be removed by the assump- 
tion that forward-angle (p,n) effective operator 
contains a spin-tensor contribution of relative 
strength 5 <~ 0.1, in addition to the Gamow- Teller 
operator, 



0(p, n ) = 

5V8^[Y 2 ({li) ® a(i)] J=1 T + (i). 



(6) 



Simple considerations of the nuclear structure 
of 71 Ga and 71 Ge suggest that the tensor oper- 
ator might be particularly troublesome for the 
71 Ga(3/2-) ^ 71 Ge(5/2-) transition to the first 
excited state. The naive description of this tran- 
sition is 



I/5/2 (n) -» 2p 3/2 (p), 



(7) 



an ^-forbidden amplitude that generates an enor- 
mous spin-tensor and vanishing Gamow- Teller 
contributions. Well-known transitions of a sim- 
ilar character, e.g. to the first excited state in 
39 K(p,n) 39 Ca, have produced discrepancies be- 
tween (p,n) and weak interaction transition prob- 
abilities of factors of ^100. The conclusion [8] is 
that it might be unwise to use the 71 Ga (p,n) 
results reliable independent measurement of 
BGT (5/2-). 

To explore this further, I did a large-basis shell 
model calculation [8] of the 71 Ga -^ 71 Ge weak 
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and (p,n) transitions. The results agree reason- 
ably with what is known experimentally: the cal- 
culated BGT (g.s.) = 0.051, compared to the ex- 
perimental value 0.087, while the calculated (p,n) 
BGT (5/2~), corresponding to the operator in 
Eq. (6), is 0.0006, in agreement with the experi- 
mental bound <0.005. However the latter result 
stemmed from a cancellation between the GT and 
spin-tensor operators comprising 0( Pin ), 



(5/2-||0 Pi „||3/2-) = 0.264- 52.23. 



(8) 



The cancellation between the second term - 
an enormous spin-tensor amplitude - and the 
GT amplitude leads to a much larger beta de- 
cay BGT(5/2~) than would be allowed by in a 
more naive interpretation of the (p,n) BGT value. 
Thus this is an explicit demonstration that pre- 
vious bounds on BGT (5/2~) are too aggressive. 

It is argued in Ref. [8] that, by relying on the 
shell model calculation of the very strong ten- 
sor amplitude in Eq. (8), a reasonable range of 
beta decay BGT (5/2 _ ) can be extracted from 
the (p,n) measurements. This results in a corre- 
sponding change in the 51 Cr cross section from 
the previous standard value 



CT ( 5i Cr) = (5.81 



to 



cr( 51 Cr) = (6.39 ±0.68) • 10 



+0.21x 
0.16) 



10- 45 cm 2 



" 45 cm 2 



(9a) 



(9b) 



where the error in Eq. (9b) represents that due 
to excited state uncertainties only. If this value 
is used in Eq. (4), one finds 

f 0.86 ±0.07 ±0.09, GALLEX . . 
\ 0.875 ±0.11 ±0.09, SAGE 1 j 

where the first uncertainty in the source exper- 
iment error while the second corresponds to the 
51 Cr cross section. Note that E ~ 1 is allowed, 
though it is certainly not demanded. It is impor- 
tant to note that the difference between (9a) and 
(9b) is one of an extended error range: all of the 
range in (9a) that is attributable to excited state 
uncertainties is allowed in (9b). It is also notable 
that the cross section uncertainty in Eq. (10) is 
comparable to the source experiment uncertainty. 
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Figure 3. The response of the pp, Be, and B 
neutrino fluxes to the indicated variations in solar 
model input parameters, displayed as a function 
of the resulting central temperature T c . (From 
Castellani et al. [13].) 



The conclusion is that the source experi- 
ments have become 7 Be neutrino cross section 
measurements. Indeed, one can express the 
GALLEX/SAGE responses to the 7 Be neutrinos 
in terms of R , independent of almost all nuclear 
physics uncertainties. This means of course that 
other tests of E^ 1 under few-atom, hot chemistry 
conditions take on added importance. Thus the 
GALLEX 71 As test discussed by Prof. Kirsten at 
this meeting is crucial. The GALLEX/SAGE re- 
sults are central to the conclusions of global anal- 
yses of the solar neutrino experiments that yield 
</>( 7 Be) < 0. 
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1.3. Solar Core Mixing of 3 He and Helio- 
seismology 

The crux of the solar neutrino problem can be 
captured in two experimental quantities. First 
the 8 B neutrino flux <j)( 8 B), which varies approxi- 
mately as T c 18 where T c is the solar core tempera- 
ture, is known to be reduced by about a factor of 
1/3 (0.47 using the new BP98 results). Naively 
this result requires a cooler sun, 

T c ~ 0.96T C SSM , 

where the superscript SSM denotes the stan- 
dard solar model result. However the flux ratio 
4>( 7 Bc)/(j)( s B), which varies as T~ w , also appears 
to be reduced relative to the standard model. 
This then requires 

T c > T C SSM 

with the extent of the increase depending on 
how strongly one wants to suppress this ratio. 
0( 7 Bc) ~ provides the best fit to the 37 C1, 
GALLEX/SAGE, and Superkamiokande results.) 
It appears that the experimental results on <j)( 8 B) 
and (j>( 7 Be)/(j)( s B) are thus in conflict, with the 
first requiring a cooler sun and the second a hot- 
ter one. 

These arguments depend on the assumption 
that neutrino fluxes will track T c as described 
above, but this appears to hold remarkably well. 
Figure 3, from Castellani et al. [13], illustrates 
this: changes due to modified nuclear cross sec- 
tions, opacity, lowered metallicity, and the solar 
age produce neutrino fluxes that track the result- 
ing T c quite accurately. This led some in the field 
to argue that no nonstandard solar model could 
produce the observed pattern of neutrino fluxes. 

Andrew Cumming and I decided to test this 
claim phenomenologically, under the assumption 
of a steady-state sun with conventional micro- 
physics producing the correct luminosity. As our 
procedures are described elsewhere [14], I'll just 
state here our basic result. There appears to 
be only one possibility for constructing a steady 
state model with a neutrino flux pattern reason- 
ably close to experiment: the solar core must 
mix on timescales of 3 He equilibration (~ few 
10 6 years) in the "elevator convection" pattern 
illustrated in Fig. 4. This mixing produces 
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Figure 4. A phenomenologically derived core con- 
vection pattern that will suppress both 4>( 8 B) 
and <^( 7 Be)/>( 8 B) [14]. The downward flow is m 
plumes, rapid and localized, requiring ~ few • 10 6 
years. This leads to out-of-equilibrium burning 
of 3 He at small r. The slow, broad, upward flow 
allows the cycle to replenish the 3 He. Typical 
upward times are ~ few • 10 7 years. 



the desired flux pattern because it modifies the 
ppl/(ppll + ppIII) and ppll/pplll branching ra- 
tios in the proper way. The former is enhanced 
because the resulting 3 He enrichment of the core 
favors the 3 He + 3 He reaction. The latter is re- 
duced because the fraction of 3 He that burns by 
3 Hc + 4 He produces 7 Be deep in the core, where 
the higher temperatures favor ppIII over ppll. 

This exercise indicated that arguments against 
nonstandard solar model solutions based on how 
neutrino fluxes scale with T c are not completely 
general. While the pattern of core mixing was 
derived phenomenologically, and not on physical 
grounds, it nevertheless has some physical ap- 
peal. The possibility of core mixing generated by 
the standard solar model overstability in the 3 He 
gradient was first discussed by Dilke and Gough. 
Roxburgh discussed a persistent convective core 
as a possible consequence of the growth of the 
SSM 3 He gradient during our sun's early convec- 
tive stage. 

Several astrophysical consequences of such mix- 
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ing were discussed, with one, helioseismology, ap- 
pearing problematic. Bahcall et al. [15] and 
Fiorcntini et al. [16] evaluated the helioseismo- 
logical consequences of replacing the SSM core 
molecular weight profile by a constant one, such 
as would occur for a continuously mixed model. 
This yielded 8% deviations in core sound speeds, 
far outside allowable bounds. My summary of 
this work is that it might be viewed as an at- 
tempt to estimate the natural scale of expected 
helioseismology changes. However, it docs not 
convincingly settle the issue because such a mod- 
ification of the SSM produces a "model" that fails 
to satisfy the equations of stellar evolution. One 
can envision that it might be more difficult to 
change the sound speed profile c(r) in a dynam- 
ically consistent model where the pressure and 
density profiles are coupled through the condition 
of hydrostatic equilibrium. 

A stronger argument, which arose from discus- 
sions at last December's ITP conference on solar 
neutrinos, is that existing helioseismology deter- 
minations of c(r) coupled with a constant molec- 
ular weight profile in the core would necessarily 
lead to an unphysical temperature profile T(r) 
that would increase away from r=0. It appears 
this conclusion is overstated: adiabatically mixed 
models lacking molecular weight gradients in the 
core generically have profiles in T(r) that are con- 
vex downward. [While this counterargument to 
the ITP discussions was originally made to me 
by Richard Epstein, similar remarks have been 
made by John Bahcall and by Doug Gough at 
this conference.] 

My view is that such mixed models, while 
quite unlikely, are not yet definitively ruled out. 
The issue of their viability is an important one, 
given the argument that a core mixed on the 3 He 
equilibration timescale is the only possibility for 
producing acceptable neutrino fluxes in a steady 
state model. Plans are underway at Los Alamos 
to evolve a series of ID models where mixing is 
included through mixing length theory, adjusting 
these in the usual way to produce the proper lu- 
minosity after 4.6 b.y. of the stellar burning. He- 
lioseismology studies performed on these models 
should then settle the issue. 



1.4. Supernova v e <-> v T Oscillations and 
Cosmologically Interesting Neutrino 

Masses 

In a Type II supernova 99% of the energy re- 
leased by the core collapse is carried off by neu- 
trinos. The initial flavor equilibrium of neutri- 
nos trapped within the core at densities p > 10 12 
g/cm 3 , coupled with the flavor-dependent decou- 
pling of neutrinos from the matter at the neu- 
trinosphere, leads to an approximate cquiparti- 
tion of energy among the flavors and to a char- 
acteristic hierarchy of temperatures. The aver- 
age energy of heavy flavor neutrinos {E VKEAVY ) ~ 
25 MeV, while {E De ) ~ 16 Mev and {E„ e ) ~ 
11 MeV. The lower values for the v e s and v e s 
reflects their stronger matter couplings due to 
charged current reactions with nucleons and to 
their greater scattering cross sections off elec- 
trons. The lower v e temperature, relative to D e , 
is due to the neutron richness of the matter near 
the neutrinospherc and resulting enhancement of 
z/ e + n^p + e~. The neutrino energy hierar- 
chy (£^ HEAVY ) > {E Ve ) > (E Ve ) appears to be 
a result independent of the details of supernova 
modelling, in contrast to the case of solar neu- 
trinos where fluxes depend on nuclear reaction 
networks. 

The spectrum of neutrinos is essentially fixed 
at the neutrinosphere, p ~ 10 12 g/cm 3 , a density 
that corresponds to a neutrino oscillation level 
crossing for Sm 2 ~ 10 4 eV 2 . Furthermore, the 
density scale height at this density will produce 
an adiabatic neutrino level crossing for sin 2 28 > 
10 . Thus the supernova neutrino spectrum pro- 
vides a unique opportunity to probe neutrino os- 
cillations. In particular, if the neutrino masses 
have a seesaw pattern 

Tfl yy o 2 2 2 

mu K Mfl m D ^ m u : m c : m t 

one can fix Mr according to the solar neutrino 
small angle solution, assuming the solar neutrino 
puzzle is due to v e — > v^. Thus u m v " ~ few 
•10~ 3 eV and "m„ T " ~ 1 eV, a value that would 
be interesting cosmologically and would induce 
v e <-» v T oscillations for supernova neutrinos. 
This is illustrated in Fig. 5. As a consequence of 
this crossing, supernova z/ e s will emerge anoma- 
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Figure 5. The three-flavor level crossing diagram 
showing two "crossings" that might be associated 
with matter enhanced oscillations of supernova 
neutrinos. 



lously hot, with an {E Ve ) ~ 25 MeV character- 
istic of heavy flavor neutrinos. An experimental 
demonstration that (E„ e ) > (E^), for example, 
would provide strong evidence for oscillations and 
possibly provide information on massive tauon 
neutrinos. 

Such an oscillation would have consequences 
for the supernova explosion, as it would enhance 
matter heating by the v e s, and on nucleosynthe- 
sis, as the increased rate for v e +n — > e~ +p would 
drive the atmosphere above the protoncutron star 
proton rich, destroying any possibility for an r- 
process. But the aspect on which I would like 
to focus is the possibility of distinctive oscillation 
signals in terrestrial supernova neutrino observa- 
tories. 

One detector of interest, despite being primar- 
ily sensitive to 9 e s, is Superkamokande. The 
usual v e signal, elastic scattering off electrons, 
will not be altered in total rate due to a v e <-> 
v T oscillation, since this rate is proportional to 
the luminosity, which is approximately indepen- 
dent of flavor. However there is some hardening 
of the spectrum of forward-scattered electrons: 
the question is whether this is enough of a sig- 



nal, given uncertainties in the supernova v e and 
^h/^h/^t/^t fluxes contributing to this forward 
scattering of electrons [17]. Perhaps more inter- 
esting is the reaction v e + w O -^ 16 F + e~, which 
produces a backward-peaked distribution of elec- 
trons that would very likely be detectable above 
the v e "background," given (E Ve ) ~ 25 MeV [18]. 
There has also been a suggestion that the 7-ray 
cascades following ^-induced spallation reactions 
on 16 might provide an attractive signal [19]. 

It seems to me that a detector of a different 
type - one flavor specific and economical to op- 
erate - might be useful in supernova watches, 
given that the characteristic time between galac- 
tic events might be ~ 35 years. An attractive 
possibility is the 1 kiloton version of the iodine de- 
tector discussed by Lande, which would be able to 
view the entire galaxy with good statistics. The 
cross section averaged over the v e flux is predicted 
to increase by a factor of ^6.6 if there is a v e <-» in- 
osculation. As the luminosity constraint leads to 
a smaller flux of (undistorted) v T s, this implies an 
increase in the iodine detector response of a factor 
of ~ 2.9, a rather dramatic signal for oscillations. 

Another possibility is the lead neutron spal- 
lation detector LAND proposed by Hargrove et 
al. [20]. The signal he discussed, single neutron 
emission, is not flavor specific nor is the cross sec- 
tion well determined. However, Fuller, McLaugh- 
lin, and I [21] recently found a rather attractive 
charge-current-specific signal in this detector, the 
emission of multiple neutrons. Due to details of 
the nuclear physics - the location of the giant res- 
onances that are strongly excited by charged and 
neutral current scattering - this channel appears 
to "filter out" neutral current events, while leav- 
ing ~ 70% of the charged current response. Thus 
if multiple neutron events are studied, LAND be- 
comes flavor specific and, due to the high thresh- 
old for reaching the giant resonances, extraordi- 
narily sensitive to the v e temperature. A v e <-> in- 
osculation would produce approximately 40 times 
the number of events that would be measured in 
the absence of oscillations. The cross section for 
detecting (E Vis ) ~ 25 MeV supernova neutrinos 
is a remarkable 4 -10 39 cm 2 . There is some work 
remaining to be done to verify these results, but 
the prospects appear quite promising. 
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This conference is a wonderful illustration of 
the power of experiments with astrophysical neu- 
trinos. Supernova neutrinos differ from atmo- 
spheric and solar neutrinos in that they impact 
the earth only once or twice a century. But they 
provide unique windows on neutrino physics, such 
as v T masses of cosmological interest. Thus it is 
probably important for us to prepare the right 
complement of neutrino observatories in antici- 
pation of the next galactic supernova. 

This work was supported in part by the U.S. 
Department of Energy. 
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